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QCD tests in r decays with optimized perturbation expansion 
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The next-to-next-to-leading order perturbative QCD corrections to R T and the higher moments of the invariant 
mass distribution in the hadronic tau decays are considered. The renormalization scheme dependence of these 
corrections is discussed. The optimized predictions are obtained, using the principle of minimal sensitivity as a 
guide to select the preferred renormalization scheme. A simplified fit is performed, using R T and R* 2 , to see how 
the use of the optimized expansion may affect the determination of the a s and the dimension six condensates 
from the experimental data. 



Recently there has been considerable interest in 
the QCD predictions for the total hadronic width 
of the r and the higher order moments of the 
invariant mass distribution in the hadronic r de- 
cays |l|-§]- These quantities have been used to 
obtain tight experimental constraints on a s and 
the condensates. In spite of the high precision 
which may be obtained in fits to the experimen- 
tal data it is important to investigate in detail 
how the results may be affected by the renormal- 
ization scheme (RS) dependence of the perturba- 
tive QCD predictions. Usually the MS scheme || 
is used to evaluate the perturbative QCD correc- 
tions. However, there is a two-parameter freedom 
in the choice of the RS in the next-to-next-to- 
leading order (NNLO), and there is no a priori 
theoretical or phenomenological reason why the 
MS scheme should be preferred. The difference 
in predictions obtained in various schemes is for- 
mally of higher order in the coupling, but numer- 
ically it may be significant, particularly at the 
energy scale of m T = 1.777 GeV. 

The QCD prediction for the R T ratio 
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where S C km = (|Kd| 2 + |K S | 2 ) w 1. The factor 
Sew = 1.0194 represents the corrections from 
electroweak interactions. The r5*°* contribution 
denotes the purely perturbative QCD correction, 
evaluated for three massless quarks. The 5^ 
contribution denotes the correction from quark 
masses (<5^ 4 « 0.009). In the case of the Rf mo- 
ments of the invariant mass distribution dT ue i/ds 
of the Cabbibo allowed decays, which are defined 
by the relation JjJ : 
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where r e is the electronic width of r, the QCD 
prediction has the form: 



R T = 3 \V u d\ 2 Sew Rq (1 + o~tt + ^svz)- 



(4) 

Rq 1 in (||) denotes the parton model prediction. 
(The S m contribution is negligible in the case of 
Rf.) The Ssvz contributions in R T and Rf de- 
note nonperturbative QCD corrections calculated 
using the SVZ approach H 



Ssvz = ^ 



Or 
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15=4,6. 



The parameters Od in Eq-@ denote vacuum ex- 
pectation values of the gauge invariant operators 
of dimension D. The cd coefficients are in princi- 
ple power series in the strong coupling constant, 
which depend on the considered quantity. 

The contribution from the D = 4 term in the 
SVZ expansion for R T may be reliably expected 
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to be small |6j], since O4 is well constrained by 
the sum rules phenomenology, and the relevant 
coefficient function starts at 0(a 2 ). However, the 
D = 6 contributions to R T is not suppressed, and 
there is little information on the value of Og. It 
was therefore proposed 0] to treat Od as free 
parameters, which are to be extracted together 
with a s from a fit to the experimental data for 
R T and the higher moments of the invariant mass 
distribution. 

The analysis reported in [Q-Q involved the Rf 
moments with (k, I) equal to (1,0), (1,1), (1,2) 
and (1,3). However, if we are interested primarily 
in the possible effect of RS dependence, we may 
simplify the discussion by considering only the 
R^ 2 moment, for which — similarly as for the R T 
— the D = 4 contribution is suppressed and there 
is significant contribution from the D = 6 term. 
Retaining in the SVZ expansion only the D = 
6 term, which appears to be a dominant source 
of the uncertainty in the nonperturbative sector, 
we obtain a simplest set of the QCD predictions 
which allows for a self-consistent extraction of a s 
and Oq from tau decays. 

The perturbative QCD corrections 5*°' and 
5p l t may be expressed as a contour integral in 
the complex energy plane with the so called 
Adler function || under the integral. (Actually 
8 P f = S$.) We have: 
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where C is a contour running clockwise from a = 
to 2 — it to = m 2 + ie away from the region 
of small \<r\. In the actual calculation we assume 
that C is a circle |er| = to 2 . The Adler function 
is defined by the relation: 



(-127r 2 )a^<V) = 3S C km[1 + 5 D . y {~a)] (7) 
r(i) 



where Ily denotes the transverse part of the vec- 
tor current correlator. The function f 12 (a/m 2 ) 
has the form: 
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The function / 00 ((t/to 2 ) has the form: 

f 00 (x) = l --x + x*- l -x\ (9) 

The NNLO renormalization group improved per- 
turbative expansion for Suy may be written in 
the form: 

6%{-v) = a(-cj)[l + na{-(T)+r2a 2 (-cT)],{lQ) 

where a = a s /ir = g 2 / [Ait 2 ) denotes the running 
coupling constant that satisfies the NNLO renor- 
malization group equation: 



da b 2 2 
cr— — --a (l + c 1 a + c 2 a ). 
da 2 



In the MS scheme we have 



„ms _ 
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1.63982 



and 1( 
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6.37101. The renormalization 



group coefficients for rif — 3 are b — 4.5, c\ = 
16/9 and cp = 3863/864 « 4.471. 

If the Adler function is expanded in terms of 
a(m 2 ), then the a dependence appears through 
the powers of ln(— a/m 2 ). The contour integra- 
tion is then straightforward and the conventional 
NNLO expansion of S^j. in terms of a(m 2 ) is eas- 



(8) 



ily obtained. However, it was observed in 11 1^] 
that one may also keep under the integral the 
renormalization group improved expression for 
the Adler function. In this case the contour inte- 
gral has to be evaluated numerically. This results 
in the essential improvement of the conventional 
expansion, corresponding to the all-order resum- 
mation of some of the corrections arising from 
analytic continuation from spacelike to timelike 
momenta. 

The QCD predictions calculated in the next- 
to-next-to-leading order (NNLO) approximation 
with massless quarks depend on two RS param- 
eters, which in principle may be arbitrary. The 
two-parameter freedom in NNLO arises because 
in each order of perturbation expansion we are 
free to choose independently the finite parts of 
the coupling constant renormalization constant. 
Different choices of the finite parts of the renor- 
malization constants result in different definitions 
of the coupling constant, which are related by fi- 
nite renormalization. (The formulas describing 
how the redefinition of the coupling affects the 



coefficients n and c 2 are coliected for exampie 
in fl3|j .) Also the dimensional QCD parameter A 
depends on the choice of the RS. In the NNLO 
there exists however a RS invariant combination 
of the expansion coefficients 
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P2 = C2 + r 2 — C\T\ 
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For the Adler function we have p 2 = 5.23783. 

The change in the expansion coefficients and 
the change in the coupling constant compensate 
each other, but of course in the finite order of 
perturbation expansion such compensation may 
only be approximate, which results in the nu- 
merical RS dependence of the perturbative pre- 
dictions. There has been intensive discussion on 
the prescriptions for making an optimal choice of 
the RS [0. One of the most attractive propos- 
tions is the choice based on the so called principle 
of minimal sensitivity (PMS) pl|, which singles 
out the scheme parameters for which the finite 
order prediction is least sensitive to the change 
of RS, similarly to what we expect from the ac- 
tual physical quantity. (It should be emphasized 
that in our case the algebraic PMS optimization 
equations |D| do not apply and a nontrivial nu- 
merical analysis of the perturbative prediction is 
required.) 

The PMS optimization has to some extent a 
heuristic character and therefore it is important 
to investigate the stability of the predictions also 
with respect to non-infinitesimal changes of the 
scheme parameters. This may be done by calcu- 
lating the variation of the predictions over a set 
of a priori acceptable schemes. A condition for 
selecting a class of acceptable schemes in NNLO 
has been proposed in 
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This condition is based on the observation, that 
the schemes with unnaturally large expansion co- 
efficients would give rise to extensive cancellations 
in the expression for the RS invariant p 2 . The 
constant I controls the degree of cancellation in p 2 
that we want to allow. In particular, for the con- 
ventional QCD expansion the PMS scheme lies 
approximately at the boundary of the region cor- 
responding to I = 2 pl| . 



A detailed discussion of the RS dependence of 
Spf has been presented in JnJ and would not be 
repeated here. An important conclusion from |l9f 
is that the contour integral resummation of higher 
order analytic continuation corrections is very im- 
portant for ensuring the RS stability of the pre- 
dictions. (This has been also discussed in p2[ . 
The instability of the conventional expansion for 



has been discussed in ]l3| , p^ , ^2] | . ) Below we 



concentrate on the RS dependence of 8, 




Figure 1. The contour plot of 5^ as a function 

of the scheme parameters r*i and c 2 , for Ai^L. = 
325 MeV. For technical reasons we use c 2 — C\r\ 
on the vertical axis instead of c 2 . The boundary 
of the region of scheme parameters satisfying the 
Eq. O is also indicated. 



The dependence of on the scheme pa- 
rameters is illustrated in figure for A^L = 

° L-l MS. — , 

325 MeV. (Similarly as in the previous work |19| 
we parametrize the freedom of choice of the RS 
by the parameters r\ and c 2 , but for technical 
reasons we use c 2 — C\T\ on the vertical axis of fig- 
ure [j] instead of c 2 .) The saddle point on figure [l] 
represents the PMS parameters. It is interesting, 
that even though the expression for 6^ has es- 



sentially non-polynomial character, the PMS pa- 
rameters are well approximated by n = and 
c 2 = t.5,02- We shall choose these values as our 
optimized parameters. (The exact PMS parame- 
ters have some dependence on the value of Ajfg. 
Also, for very large values of A-p^ the pattern 
of RS dependence is more complicated than this 
shown in figure |[ However, for all values of A-p-j 
the RS dependence in the vicinity of r% = and 
C2 = 1.5,02 is very small.) The set of scheme pa- 
rameters, that involve the same — or smaller — 
degree of cancellation in p 2 than our preferred pa- 
rameters, satisfies the condition ( p"3| ) with I = 2. 
The boundary of this set is indicated on the fig- 
ure [l| By calculating the variation of the predic- 
tions over this set of scheme parameters we may 
estimate in a quantitative way the sensitivity of 
the NNLO predictions to the change of RS. Thus 
obtained estimate may then be compared with 
a similar estimate for other quantities, for which 
the NNLO predictions are known. It should be 
noted that although the MS parameters lie out- 
side of the "allowed" region shown in the figure [l], 
the numerical value of the prediction in the MS 
scheme is close to the lowest value attained in the 
"allowed" region. 

It is interesting to note, that the RS depen- 
dence pattern of is quite similar to the pattern 
found in [|l9| for 8 °^, despite the fact that the con- 
ventional expansions for these quantities appear 
to be quite different. Indeed, in the conventional 
expansion for S^f we have r^ IS = 5.2023 and 
r^S = 26.3659, which gives p 2 = -5.4757, while 
for 8$ we have rf s = 3.5795 and rf^ = 4.3441, 
which gives pi = —10.3614. This seems to sug- 
gest that the improved predictions obtained with 
contour integral expressions are more natural, re- 
flecting the common origin of the two corrections 
in a better way. 

In figure || we show the RS dependence of the 
next-to-leading order (NLO) prediction for S^f. 
By numerical optimization we find that in NLO 
r PMS w _Q.64. (Again, this depends to some ex- 
tent on A-p-j, but this dependence has negligible 
effect on the numerical value of the prediction.) 

In figure || we show the NNLO PMS predic- 
tions for 81+ as a function of m r /Al^_, together 
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Figure 2. NLO prediction for 8^ as a function of 
the scheme parameter r\, for Ai-L = 325 MeV. 



with the minimal and maximal values obtained 
by varying the scheme parameters within the re- 
gion determined by the condition ( |l3| ) with 1 = 2. 
We see that the NNLO predictions for 5*1, ob- 
tained by numerically evaluating the countour in- 
tegral expression (Q), are free from potentially 
dangerous RS instabilities, even for large values 

(3) 

of A-^— . This situation is similar to that cncoun- 

MS 

tered for o"*°*. For comparison we also show the 
PMS predictions obtained in the next-to-leading 
order (NLO). We see that RS dependence of the 
NNLO expression within the region defined by 
the condition (|l^) is smaller than the difference 
between NNLO and NLO PMS predictions. In 
a separate figure we show the NNLO and NLO 
predictions in the MS scheme (figure ||) . 

In order to see, how the optimization of the 
scheme parameters affects the fits to the exper- 
imental data, we first test the accuracy of the 
approximation in which one only retains the Oq 
contribution in the SVZ expansion. To this end 
we make a fit of a s and Oe in the MS scheme 
and compare the results with the fit performed 
by ALEPH ||, in which the 4 , O e and O s con- 
tributions have been taken into account in the 
(1,0), (1,1), (1,2) and (1,3) moments. To make 




Figure 3. The optimized predictions for 8l£, as a Figure 4. as a function of m T / A^—L, obtained 



function of m T /Aj-L-, obtained in NNLO (upper 
solid curve) and NLO (lower solid curve). The 
dashed lines indicate variation of the predictions 
when the scheme parameters are changed within 



the region satisfying the Eq. [13 



the fits we use the following expressions: 

R T = 3 x 1.0194 (0.991 + <5* ot - 3.75 6 ), (14) 



and 
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210 (l + ^ ot -3.75 6 ) 

If we take, following ALEPH §, R T = 3.645 ± 
0.024 and D\ 2 = 0.0570 ± 0.0013, we obtain from 
the fit in the IIS scheme a s {M%) = 0.1209 ± 
0.0013 and O e = -0.0010 ± 0.0012. This appears 
to be remarkably close to the values 0.121 and 
-0.0016 obtained in the full fit by ALEPH. This 
gives us confidence that the u O@ approximation" 
captures the essential features of QCD corrections 
in t decays. 

Let us now study how the RS dependence may 
affect the fit to the experimental data. Let us 



P* " "MS 7 

in the MS scheme in NNLO (upper curve) and 
NLO (lower curve). 



take R T = 3.635±0.016, which is a weighted aver- 
age of three possible determinations [^0| , involv- 
ing B e = 0.1783 ± 0.0008, = 0.1735 ± 0.0010 
and r r = (291.0 ± 1.5) x 10~ 15 sec. (These are the 
so called "our fit" values. Using a set of "our av- 
erage" values given in [EO we would obtain R T = 



3.643.) Let us also take 



D]r = 0.0559±0.0007, 



which is the most precise published value up to 
date. Using the NNLO PMS expression we ob- 
tain then a a (M§) = 0.1188 ± 0.0008 (a s (m 2 T ) = 
0.330 ± 0.008) and O e = -0.0021 ± 0.0006. Per- 
forming the same fit, but using now the NNLO 
MS expression, we obtain a a {M%) = 0.1198 and 
6 = -0.0020. We see that the value of the 
condensate practically does not change, but the 
value of the strong coupling constant is reduced 
by amount approximately equivalent to the ex- 
perimental error. 

It is also of some interest to repeat the fit us- 
ing the NLO PMS expression. We then obtain 
a s (-M§) = 0.1210 and the same value for Oq as 
in NNLO. Taking the difference of the NNLO and 



NLO PMS fits is perhaps the best way of estimat- 
ing the accuracy of the perturbative contribution 
in this problem, since we found the QCD correc- 
tions to be quite stable with respect to change 
of the scheme. We see that thus obtained un- 
certainty of the fitted value of a s is of the order 

0. 0022. i. e. it is quite large, compared for example 
to the experimental uncertainty. (This estimate 
of uncertainty involves only the perturbative un- 
certainty — to have estimate of theoretical un- 
certainty for the total QCD prediction one should 
also consider the accuracy of the SVZ expansion 
itself.) It should be noted that the NNLO-NLO 
difference is strongly RS dependent so it is essen- 
tial to optimize the choice of the scheme before 
comparing the predictions in successive orders. 

Concluding, the renormalization scheme de- 
pendence of the perturbative QCD corrections 
to R T and to the R^. 2 moment of the invariant 
mass distribution in hadronic tau decays has been 
studied in detail. The optimized predictions have 
been obtained using the principle of minimal sen- 
sitivity as a guide to select the preferred renor- 
malization scheme. The stability of the predic- 
tions obtained via the contour integral expres- 
sion (^) has been verified, using a specific con- 
dition to eliminate the schemes that have unnat- 
urally large expansion coefficients. However, the 
difference between predictions in the convention- 
ally used MS scheme and the optimized predic- 
tions obtained in the PMS scheme was found to 
be phenomenologically significant. Also, the dif- 
ference between the NNLO and NLO predictions 
in the PMS scheme was found to be significant, in- 
dicating perhaps that the uncertainty in the per- 
turbative expression is larger than previously ex- 
pected. 

REFERENCES 

1. F. LeDiberder and A. Pich, Phys. Lett. B 289 
(1992) 165. 

2. ALEPH Collab., D. Buskulic et al, Phys. Lett. 
B 307 (1993) 209. 

3. P. Reeves, in Proceedings of the XXXth Ren- 
contres de Moriond Conference "'95 QCD 
and High Energy Hadronic Interactions," Les 
Arcs, Savoie, France, 1995, edited by J. 



Tran Thanh Van (Editions Frontieres, Gif- 
sur-Yvette, 1995), p.235. 

4. CLEO Collab., T. Coan et al, Phys. Lett. 
B 356 (1995) 580. 

5. W.A. Bardeen, A.J. Buras, D.W. Duke and 
T. Muta, Phys. Rev. D 18 (1978) 3998. 

6. E. Braaten, Phys. Rev. Lett. 60 (1988) 1606, 
ibid. 63 (1989) 577, S. Narison and A. Pich, 
Phys. Lett. B 211 (1988) 183, E. Braaten, 
S. Narison and A. Pich, Nucl. Phys. B 373 
(1992) 581. 

7. M.A. Shifman, A.I. Vainshtein and V.I. 
Zakharov, Nucl. Phys. B 147 (1979) 385, 448, 
519. 

8. C.S. Lam and T.M. Yan, Phys. Rev. D 16 
(1977) 703, K. Schilcher and M.D. Tran, 
Phys. Rev. D 29 (1984) 570. 

9. S.L. Adler, Phys. Rev. D 10 (1974) 3714. 

10. S.G. Gorishny, A.L. Kataev and S.A. Larin, 
Phys. Lett. B 259 (1991) 144, M.A. Samuel 
and L.R. Surguladze, Phys. Rev. D 44 (1991) 
1602. 

11. A. A. Pivovarov, Z. Phys. C 53 (1992) 461. 

12. F. LeDiberder and A. Pich, Phys. Lett. B 286 
(1992) 147. 

13. P.A. Raczka, Phys. Rev. D 46 (1992) R3699. 

14. P.A. Raczka, in Proceedings of the XV In- 
ternational Warsaw Meeting on Elementary 
Particle Physics, Kazimierz, Poland, May 
1992, edited by Z. Ajduk, S. Pokorski and 
A.K. Wroblewski, World Scientific, Singa- 
pore, 1993, p. 496. 

15. P.M. Stevenson, Phys. Rev. D 23 (1981) 2916. 

16. A. Dhar, Phys. Lett. 128B (1983) 407. 

17. For the summary of early contributions see 
D.W. Duke and R.G. Roberts, Phys. Rep. 120 
(1985) 275. 

18. P.A. Raczka, Z. Phys. C 65 (1995) 481. 

19. P.A. Raczka and A. Szymacha, Z. Phys. C 70 
(1996) 125. 

20. Review of Particle Properties, Particle Data 
Group, L. Montanet et al, off-year partial up- 
date (January 1996). 



